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Abstract 

We consider weighted Reed-Muller codes over point ensemble Si x 
• ■ • x S m where Si needs not be of the same size as Sj. For m = 2 we 
determine optimal weights and analyze in detail what is the impact of 
the ratio |Si|/|S2| on the minimum distance. In conclusion the weighted 
Reed-Muller code construction is much better than its reputation. For a 
class of affine variety codes that contains the weighted Reed-Muller codes 
we then present two list decoding algorithms. With a small modification 
one of these algorithms is able to correct up to 31 errors of the [49, 11, 28] 
Joyner code. 

Keywords. Affine variety codes, list decoding, weighted Reed-Muller 
codes 

1 Introduction 

Weighted Reed-Muller codes were introduced by S0rensen in |28j . In his paper 
he demonstrates that they are subcodes of q-axy Reed-Muller codes of the same 
minimum distance and it is therefore not surprising that not much attention has 
been given to them since. In the present paper we consider the above two code 
constructions in a slightly more general setting as we allow any point ensemble 
S = Si x • • • x S m , Si,...,S m C F q . Other authors have considered g-ary 
Reed-Muller codes in this setting, but nobody seems to have recognized that 
for such point ensembles weighted Reed-Muller codes are often superior. We 
shall derive a number of results regarding their efficiency and define what we 
call optimal weighted Reed-Muller codes in two variables. 

We argue that the dual codes are exactly as efficient and that they can be de- 
coded up to half the designed minimum distance by known decoding algorithms. 
We then turn to the decoding of weighted Reed-Muller codes. The first decoding 
algorithm that we present utilizes the fact that the codes under consideration 
can be viewed as subfield subcodes of certain Reed-Solomon codes. This algo- 
rithm is a straightforward generalization of Pellikaan and Wu's list decoding 
algorithm [21]. The second decoding algorithm that we present is a more di- 
rect interpretation of the Guruswami-Sudan list decoding algorithm. We are by 
no means the first authors to consider such an approach for multivariate codes 
(see [21], PQ, [2J). Our contribution is that we develop a method for deriving im- 
proved information on how many zeros of prescribed multiplicity a multivariate 



polynomial can have given information about its leading monomial with respect 
to the lexicographic ordering. Using such information and allowing the decod- 
ing algorithm to perform a preparation step we develop an improved algorithm. 
For some optimal weighted Reed-Muller codes the first decoding algorithm of 
the paper is quite good, for others the latter is the best. 

Weighted Reed-Muller codes are examples of a particular class of affine variety 
codes. Whenever possible we state our findings for this more general class of 
codes. As a bonus we find that when equipped with a small trick the subfield 
subcode decoding algorithm can decode the Joyner codes [HJ Ex. 3.9] beyond 
its minimum distance even though till now this code has resisted even minimum 
distance decoding. 

2 A class of affine variety codes 

Given 

S = S 1 x---xS m = {P 1 ,...,P lsl } 
write n = \S\ and consider the evaluation map 

ev 5 : F q [X lt ...,X m ]-> FJ, ev 5 (P) = (F(J\), . . . , F(P n )). 

Let 

\0<ij < 1^-1^ = 1,...,™} 

and define the affine variety code 

E(M,S) = Span Fa {ev s (M) | M e M}. 

Throughout the paper we use the notation s, = \S%\ for i = 1, . . . , to. If not 
explicitly stated we shall always assume that the enumeration is made such 
that si > • • • > s m holds. In the special case that S\ — ■ ■ ■ — S m we write 
S = S x • • • x S and s — \S\. We first show how to find the dimension of the 
code. 

Proposition 1. The dimension of E(M,S) equals |M|. 

Proof. We only need to show that 

{evspqv . .,X%) \ 0<ij< Sj,j = l,...,m} 

constitutes a basis for F™ as a vectorspace over F q . For this purpose it is 
sufficient to show that the restriction of ev^ to 

{G(X 1 ,...,X m )\deg Xi G<s h i=l,...,m} (1) 

is surjective. Given (ai, . . . , a n ) € F™ let 

n m / Y \ 

=5>n n \r\ 

Here, we have used the notation P v = (P^ , . . . , P^), V = 1, ...,n. It is 
clear that eys(P) = (01, ... ,a„) and therefore evs : F q [Xi, . . . ,X m ] — > F™ is 
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surjective. Consider an arbitrary monomial ordering. Let R(Xi, . . . , X m ) be 
the remainder of F(X\, . . . , X m ) after division with 

{ JJ (X t - a),... , J] (X m -a)}. 

Clearly, F{Pi) — R{Pi) = flj, i = 1, . . . , n. Hence, the restriction of evs to ([!]) 
is indeed surjective. □ 

We next show how to estimate the minimum distance of E(M,S). The 
Schwartz-Zippel bound [251 1301 IS] is as follows: 

Theorem 2. Given a lexicographic ordering let the leading monomial of F(X\, . . . , X m ) 
be Xl 1 ■ ■ ■ X^ . The number of elements in S = Si x • • • x S m that are zeros of 
F is at most equal to 

HS2 ■■■Srn+ Sl«2S3 ' ' ' s m H h Si ■ ■ ■ s m _ii m . 

The proof of this result is purely combinatorial. Using the inclusion-exclusion 
principle it can actually be strengthened to the following result which is a special 
case of the footprint bound from Grobner basis theory: 

Theorem 3. Given a lexicographic ordering let the leading monomial of F(X\, . . . , X m ) 
be Xl 1 ■ ■ ■ X'^ 1 . The number of elements in S = Si x • • • x S m that are zeros of 
F is at most equal to 

n~{si- i 1 )(s 2 - i 2 ) ■ ■ ■ (s m ~ i m ). 

Proposition 4. The minimum distance of E(W.,S) is at least 

min{(si - ii){s 2 -i 2 ) • ■■ (s m - im)\X^ ■ ■ ■ X 1 ™ e M}. 

The bound is sharp if for every M el all divisors of M also belong to M. 

Proof. The first part follows from Theorem [3j To see the last part write for 
i = 1, . . . ,m, Si = {bi \ . . . , The polynomial 

F{Xi,...X m ) = \\X{{X v -bf) 

v=\j = \ 

has leading monomial X 1 ^ ■ ■ ■ X^ with respect to any monomial ordering and 
evaluates to zero in exactly n — (si — ii)(s 2 — i 2 ) • •• (s m — i m ) points from S. 
Finally, any monomial that occurs in the support of F is a factor of X} 1 • • • X]^ . 

□ 

3 Weighted Reed-Muller codes 

The first example of codes E^M, S) that comes to mind are the g-ary Reed- 
Muller codes RM 5 (m, m). They are defined by choosing 

Si = • • ■ = S m = Fq, 
M = {X{ 1 ---X^ \ii + ... + im <u}. (2) 
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S0rensen in [28] modified the above construction by instead letting 

M = {XY ■ ■ ■ X% I toiti + • • ■ + w m i m < u} (3) 

where w% , . . . , w m are fixed positive numbers. The resulting codes are called 
weighted Reed-Muller codes. In the same paper S0rensen argues that there is 
actually no point in considering ^ rather than ^ as every weighted Reed- 
Muller code is contained in a code RM,(ii,m) which has the same minimum 
distance. In the present paper we allow Si, ... , S m to be any subsets of F q . As 
we shall demonstrate, in such a general setting replacing ^ with ^ may result 
in much better codes. In other words, the concept of weighted Reed-Muller 
codes actually makes a lot of sense. We start with a motivating example. 

Example 5. In this example we construct codes over F 16 of length n = 64. 
First let S — Si x S 2 be such that s± = s 2 = 8. Define, 

M = {Xl'Xi? I < h, (3 < 7, it + i 2 < 7}. 

The code E(M.,S) is of dimension 36 and minimum distance 8. Letting instead 
S = Si x S 2 where \S±\ = 16 and \S%\ =4 we consider the following two sets of 
monomials 

M = {XI 1 X? I < h < 15, < i 2 < 3, h + i 2 < 11}, 
M" = {X^Xl 2 I < ii < 15, < i 2 < 3, h + 2i 2 < 14}. 

The code _E(M', S) is of dimension 42 and minimum distance 8 whereas the code 
E(Wl / ' 1 S) is of dimension 48 and minimum distance 8. 

The above example illustrates two facts. Firstly, choosing the Si's to be of 
different sizes may be an advantage. Secondly, using a weighted degree rather 
than the total degree when choosing monomials may result in better codes. It 
is time for a definition. 

Definition 6. Let Si, . . . , S m C V q and consider positive numbers w\, . . . , w m , u. 
Let 

M = {XI 1 ■ ■ ■ X^ I < it < s t — l,t = 1, . . . , rn, and W\ii H + w m i m < u}. 

The corresponding code E^M, S) is called a weighted Reed-Muller code and we 
denote it by RM(S%, . . . , S m , u, uii, . . . , w m ). As is often done we shall refer to 
weighted Reed-Muller codes with Si = ■ ■ ■ = S m and w\ = ■ ■ ■ = w m as q-ary 
Reed-Muller codes. 

We start by taking a closer look at the case of two variables. According to 
Theorem [4] the minimum distance of RM(Si, S 2 , u, u>i, w 2 ) equals 

min{(si-ii)(s 2 -i 2 ) | H,«2 € N, < ii < Si-1,0 < i 2 < s 2 -l,wiii+w 2 i 2 < u} 
> min{(si - ii)(s 2 - i 2 ) \ ii,i 2 e Q, < i\ < si - 1, 

< i 2 < s 2 - 1, wiii + w 2 i 2 = u}. (4) 

Substituting i 2 — (u — wiii)/w 2 into (si — ii)(s 2 —i 2 ) we get a concave function 
(a parabola). Hence, the minimal value of (si — ii){s 2 — i 2 ) under the condition 
in Q is either attained for ii as small as possible or for ii as large as possible. 
Given a weight wi and a positive number u we seek w 2 such that (si — ii)(s 2 — i 2 ) 
is the same for ii minimal and maximal under the condition in El). 
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Figure 1: The situation in the proof of Proposition [7] 



Proposition 7. Let S2 < Si 6e positive integers. Given fixed positive numbers 
W\ and u assume W2 is chosen to be the positive number such that (s\— ii)(s2— 12) 
attains the same value whenever i\ is minimal or is maximal under the condition 

w x ix + w 2 «2 = u, 
0<«i<si-l, 0<i 2 <s 2 -l. 

We have 

if < u < ( Sl - %)wi 
if Oi - f^Vi < u < (si - l)wi 
£/ (si - l)wi < u < (si - l)wi + (s 2 - 1)^2- 

(5) 

Proof. The proposition is illustrated in Figure [T] for the case of Si = 18 and 
s 2 = 6. 

We concentrate on the situation where 

(si < it < (si — l)u?i 

and leave the other two simpler cases for the reader. Write u/w\ = s% — 5 with 
s i/ s 2 > 8 > 1. The maximal value of i\ is u/wi in which case i 2 = 0. So 
for ii maximal (si — «i)(s2 — ^2) = ^^2- We seek i\ minimal such that with 
i2 = S2 — 1 we get (si — ii)(s 2 — i 2 ) = <5s2- We find i\ = si — <5s2 which is indeed 
a non-negative number. Hence, u?2 must satisfy 

wi(s\ — <5s 2 ) + w 2 {s2 - 1) = — (5) 

Wl 1 

U)2 5 

W2 W\S\ — U 

□ 

Proposition [7] justifies the following definition. 

Definition8. If s\, S2, u, w\, W2 satisfy |5|) then the code RM(S\, S 2 ,u,w\,W2) 
is called an optimal weighted Reed-Muller code (in two variables). 




S2/S1 
Wi/{w 1 s 1 - u) 
1 
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The next proposition estimates the minimum distance of any weighted Reed- 
Muller code RM(Si, S2, u, Wi,w 2 ) (optimal or not). 

Proposition 9. Consider RM(S\, S2,u,wi,w 2 ) with s 2 < s\. Write p = w\jwi 
and let d be the minimum distance. 
Ifp<^ then 

d>s 2 (si ), ifu < (si - (6) 

Wl 

d>s 2 - ~ ~ ^ ~ — , if (si - l)w x <u<( Sl - l) Wl + (s 2 - l)w 2 . (7) 
w 2 

Iff<p<\ then 
u 

d>(s 2 )sx, ifu < (s 2 - l)w 2 , (8) 

w 2 

u — (s 2 — IW2 „, 1 N 

d> Sl ^ >—L if{s2-l)w2<u<{s 1 --)w 1 , 9 

Wi p 

u 1 

d>(si )s 2 , if (si )wi < U < (si - l)u;i, (10) 

wi p 

d>s 2 - U - ( ' Sl - — , if (sx - l) Wl <u<( Sl - l) Wl + (s 2 - l)w 2 . (11) 
w 2 

If 1 < p then 

u 

d>(s 2 )si, if u < (s 2 - l)w 2 , (12) 

w 2 

d>s x - U - (§2 - 1)W2 , if (s 2 - l)w 2 <«<(«i- l)u>i + (* 2 - l)w 2 - (13) 

Equality holds in Q) 7 Q), arid |'iS), respectively, if the expression is an 

integer. Equality holds in (T6t) and (10) if u/w\ is an integer. Finally, equality 
holds in and if ujw 2 is an integer. 

Proof. The task is to determine under the various conditions of the proposition 
whether (si — ii)(si — i 2 ) is minimized for i\ minimal or maximal. The corre- 
sponding values of i\ and i 2 are then plugged in to give Q,- • • ,( 13 ). To find out 
if i\ should be chosen minimal or maximal we use the information from Propo- 
sition If p < s 2 / Sl the minimum is always attained for i-y maximal. If 1 < p 
then the minimum is always attained for i\ minimal. In the case s 2 j s\ < p < 1 
the minimal is attained for i\ minimal when u < u' and is attained for i\ max- 
imal when u > u' . Here, u' is a number that we determine below. It is clear 
that 

{si -)wi < u' < (s! — l)wi 

s 2 

and therefore u' is the number such that 

v! u' - (s 2 - l)w 2 
(si )s 2 = si . 

Wl Wl 

Solving for v! gives 

( X \ 

U = W1S1 - w 2 = (Si )Wi. 

P 

□ 



G 



Proposition [9] also allows us to state general bounds for the minimum dis- 
tance of RM(S'i, . . . , S m , u, wi, . . . , w m ) in two important cases. Observe, that 
in particular the following proposition can be applied when Wi — Yii^j s j ■ 

Proposition 10. Assume s\ > • • • > s m , and let u be a number < u < 
(si - l)wi H h (s m - l)w m . If 

w i < ^ <...< ™" (14) 



holds then write 

u = (si - l)w\ H h (sj-i - l)iy t _i + a t w t 

where < a t < s t ~ 1. TTie minimum distance of RM(Si, . . . , S m ,u,wi, . . . ,w m ) 
satisfies 

m 

d>(s t - at) } J s 4 

i=i+l 

with equality if at is an integer. 
If Wi > ■ ■ ■ > w m then write 

u = (s m - l)w m H h (s 4 _i - + a t w t 

where < a t < s t — 1 . TTie minimum distance of RM{S\ , . . . , S m , u, wi, . . . , w m ) 

> (s t - at) Y[ s i 

i=l 

wit/i equality if at is an integer. 



Proof. We only prove the first part. Assume (14) holds. Let ix, . . . , i m £ Q be 



chosen such that (si — «i) • • • (s m — i m ) is minimal under the conditions 

Will H h W,„?', Tl = U, 

< i\ < si - 1, . . . , < i m < s m - 1. 

For integers c,d with l<c<d<77iwe have w c /wd < Sd/s c . Note from 
Proposition [7] that Sd/s c is the smallest possible ratio of w' c /w' d for an optimal 
code RM(S' C , Sd, i c w c +idWd, w' c , w d ). Therefore, under the condition that i c w c + 
idWd is fixed and < i c < s c — 1, < id < Srf — 1 the minimal value of (s c — 
ic){sd — id) is attained for id minimal. The result now follows by induction. □ 

In the remaining part of this section we restrict solely to the case of two 
variables. As shall be demonstrated in this situation almost all weighted Reed- 
Muller codes outperform the corresponding g-ary Reed-Muller codes. Before 
getting to the analysis let us consider an example. 

Example 11. Consider optimal weighted Reed-Muller codes RM(Si, S2, u, vi\, W2) 
(see Definition^. Choosing (si,S2) from the set 

{(32, 32), (64, 16), (128, 8), (256, 4), (512, 2)} 
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d 
n 



32 x 32 
64 x 16 
128 x 8 
256 x 4 
512 x 2 




Figure 2: Performance of the codes in Example [XT] 



gives five different classes of codes all of length n = 1024. Observe that the 
first class of codes is similar to q-ary Reed-Muller codes as the optimal choice 
of W\,W2 is Wi = W2 whenever si = S2- The codes are defined whenever the 
field under consideration contains at least si elements. Hence, the first class of 
codes is defined over any field F q with q > 32, the second class over any field 
F q with q > 64, the last class of codes over any field F q with q > 512. In 
particular all classes of codes are defined over -F512. In Figure [77] we compare 
their performance. It is clear that the second class of codes outperforms the first 
class for higher dimensions, whereas the last three classes of codes outperform 
the first class for any dimension. 

Below we investigate in detail how well general optimal weighted Reed- 
Muller codes RM(5i, 52, u, Wi, W2) perform in comparison with q-aiy Reed- 



Muller codes RM(5, 5, it', 1, 1). Here, we assume that S1S2 



Recall, 



from Proposition [7] that the description of the weights used in the optimal 
weighted Reed-Muller codes involves three cases depending on the value of u. 
Choosing in the following without loss of generality w\ = 1 we shall refer to 
u < s\ — (S1/S2) as region I, s± — {S1/S2) < u < s± — 1 as region II, and finally 
Sx — 1 < u < (si — 1) + w 2 {s2 — 1) as region III. Proposition 12 Proposition 13 



and Proposition |14| respectively, takes care of region I, region II, and region 
III, respectively. In Proposition [12] we will to ease the analysis make the small 
restriction that S2 \ si and that s\ \ us2- Furthermore, in all three propositions 
we assume that u is an integer. We stress that when such assumptions do not 
hold then the formulas to be presented are still very close to be true. What 
we will learn is that the codes RM(5i, 52, u, wi = 1, u> 2 ) always outperform the 
codes RM(5, 5, u' , 1, 1) provided that si > 4s2- Furthermore, for si — s/s 2 < u 
such a result holds in the general situation s± > S2- 

Proposition 12. Consider integers Si, S2 with 1 < S2 < s±. Let u be an integer 
with u < si— {sij s%). Assume s\/ S2 and US2/ s\ are integers and that S1S2 = s 2 
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for some integer s. Let W\ = 1 and W2 = S1/S2 (that is, w\ and w% are chosen 
as in Proposition^. The code RM(S\, S2, u, wi, W2) is of dimension 

- (u 2 — +u ) + u— + 1 
2 si si 

and any RM(S, S, w', 1, 1) of the same or larger minimum distance is of dimen- 
sion at most 

I(£V + 3u # + 2). 
2 si y si 

For si > 4s2 t/ie code RM(S\, S2, u, Wx, W2) is the better one. 
Proof. The dimension of RM(5i, S2, u, ui\,W2) is 

us 2 /si . 
^ I US 2 
> —S 2 + + 1 

z— 1 

and the minimum distance is S2(si — u). Assuming v! = U \J~^ is an integer, 

the code RM(>!?, S, u' , 1, 1) is of minimum distance S2(si — u). This code is of 
dimension 

2 si V s i 

□ 

Proposition 13. Consider integers s\ and S2 with S2 < s±. Let u be an integer 
with 

S\ < U < Sl — 1. 

S2 

Assume Sis 2 — s 2 for some integer s. Let w\ = 1 and W2 = S\—u (that is, w\ 
and W2 are chosen as in Proposition^. The dimension of RM(S\, S2, u, Wi, W2) 
equals 

giCgl - u ) , s 2 (si -u) 
sis 2 + s 2 2 ■ ( 15 - ) 

Lf u> s\ — s/s2 then any code RM(S, S, v! , 1, 1) 0/ £/ie same or larger minimum 
distance is of dimension at most 

(s x - u)s 2 ((s! - u)s 2 - 1) 
S1S2 „ 



which is less than (15) for s 2 < s\. Lfu< s±— s/s 2 then any code RM{S, S,u" ,1,1) 
of the same or larger minimum distance is of dimension at most 

l(^ + 2)(^ + l). 



TTiis number is smaller than the value of (15) for S\ > 4s2 and equal if si = 4s2. 

Proof. Consider the first code which is of minimum distance (si — u)s2- For 
?2 = S2 — 1 the value i\ such that w\i\ + w 2 i 2 = it is i\ = u — (sa ~ l)( s i — u )- 
Therefore the dimension equals 

S2(ll + 1) + 22 ( S l - W ) Z = S 1 S 2 - + S2 - 

i=l 
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If (si — u)s2 < s u > S\ — s/s 2 then for v! = 2s — 1 — sis 2 + US2 the code 
RM(5, S, vf, 1, 1) is of the same minimum distance. This code is of dimension 

2 ■ (si - tt)s 2 ((si - n)s 2 - 1) 
s - 2^ « = sis 2 ^ ■ 

The dimension of the first code exceed the dimension of the latter code by 
(si — u — l)fs 2 (si — — ^2) which is a positive number for u < s\ — 1 and 
equals zero for u = S\ — 1. 

If (si — m)s2 > s w < si — s/s2 then imagining that s divides S2{s\ — u) the 
code RM(S, 5, it", 1,1) with 

(s - u")s = s 2 (s! - u) u" = — 

s 

is of the same minimum distance as RM(Si, S2, u, wi, W2). The dimension equals 

u" + l 



E 

i=l 



1 / S 2 U „\ / S 2 U 



Subtracting this expression from ( 15 ) one gets a concave function (a parabola) 
in u. Therefore the smallest value of the difference is attained either for u = 
s i ~ s i/ s 2 or f° r u = Si — s/,s 2 . Plugging in the first value and substituting 
si = xs, S2 — s/x one finds that the resulting function is zero for x = 2 and 
positive for x g]2; s[. Plugging in the latter value is not needed as we already 
know from the first part of the theorem that here the difference is positive. □ 

Proposition 14. Consider integers si and S2 with 1 < S2 < si and sis 2 = s 
where s is an integer. Let u be an integer with s\ — 1 < u < (si — 1) + (s 2 — 1). 
Let wi = W2 = 1 (that is, w\ and W2 are chosen as in Proposition^^. There 
is a Reed-Muller code over S x S of the same minimum distance and the same 
dimension. 

Example 15. This is a continuation of Example \11\ Consider the graph in 
Figure [7I| First we take a look at the optimal weighted Reed-Muller codes cor- 



responding to (si,s 2 ) = (64,16). For these codes region I (Proposition 12) 
corresponds to rates k/n below approximately 0.5. As 64 = 4 • 16 we expect the 
optimal weighted Reed-Muller codes to behave very much like the correspond- 
ing Reed-Muller codes in this region, which is indeed what the graph reveals. 
Considering values o/(si,S2) with si > 4s 2 , when S1/S2 increases the rates cor- 
responding to region I defines a smaller and smaller interval (starting of course 
still with rate equal to 0). The improvements in region I increases, but a more 
important contribution for the codes to become better and better is that region 
11 takes over at smaller rates. Similarly, the interval of rates corresponding to 



region III (Proposition 14) becomes smaller and smaller. This is the interval 
where the optimal weighted Reed-Muller codes are (again) as bad as the q-ary 
Reed-Muller codes. For (si,s 2 ) = (64,16) this last mentioned interval starts 
at approximately 0.875. Already for (si, S2) = (128, 8) the starting point of the 
interval is around 0.97. 

Proposition [l2j Proposition [13] and Proposition [14] tell us that whenever 
si > 4s2 then the optimal weighted Reed-Muller codes outperform the Reed- 



Muller codes coming from S x S. Example 11 further suggests that from that 
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point further increasing s% and decreasing s 2 can only help. The following three 
propositions together confirm this observation. As in previous propositions we 
will need to make a few assumptions on the codes that we consider. Again we 
stress that when such assumptions do not hold then the formulas to be presented 
are still very close to be true. 

Proposition 16. Consider positive integers S\,s 2 , s^, s' 2 with S1S2 = s'iS 2) s[ > 
Si > S2 > s' 2 and s\/ s 2 > 2. Let u, < u < s% — Si/s 2 be an integer with S\\us2 
and consider the optimal weighted Reed-Muller code 

RM(S 1 ,S 2 ,u,w 1 = l,w 2 = 81/82). (16) 

If for an integer v! with s[\u' ' s' 2 RM(S{, S 2 , vf, w[, w 2 ) is an optimal weighted 



Reed-Muller code of the same or smaller minimum distance as that of (16) then 
the latter code is of dimension at least that of the first one. 



Proof. If the latter code is of minimum distance close to that of ( 16 ) then it 
belongs to region I or II. The codes are of minimum distance d = s 2 (si — u) and 
d! = s' 2 (s' 1 — u'), respectively. Hence, u' > us 2 /s' 2 . If the latter code is in region 
I the improvement in dimension is at least 

_ ,1 £2 s 2 s 2 
U( 2 + Sl 2s> 2 s'/ 

which is positive when s 2 + si/2 < s^ + s^/2. Writing s[ = fisi this corresponds 
to 

^ 2 (si/2) + fj,(-s 2 - Si/2) + s 2 > 0. (17) 
The left side is a convex parabola with roots /1 = 1 and = 2s 2 /s\. The as- 



sumption Si/s 2 > 2 therefore guarantees that (17) holds for all \x > 1. 
Assume next that the latter code in the proposition is in region II. The improve- 
ment can be calculated to be 

« 2 (-^) + ^ + f-^) + (^-^ + 4-D (18) 

which is a concave function in u. Our assumptions give s± — s\/s 2 < u < 
S\ — S\/s 2 and therefore it is enough to plug u — si — Si/s' 2 and u — si — Si/s 



into ( 18 ) and then to check that the resulting values are positive. The first value 



is positive if 



+ Af2(1 + 2T 2 + T } + ~ I " S2) + 82 (19) 



is positive. The roots of this function in \i are 0, 2s 2 /si, and S2. Hence, (19) is 



indeed positive for e]l,s 2 [. When u = — si/s 2 is plugged into (18) we get 



S 2 SiS 2 SiS 2 

S2 

ix 2 2n 

which is positive for ji > 1. □ 

Proposition 17. Consider positive integers s\,s 2 , s' 1} s 2 with s±S2 = s'iS 2) s[ > 
si > s 2 > s' 2 , and s\/s 2 > 2. Let u, s\ — s±/s 2 < u < si — 1 be an integer and 
consider the optimal weighted Reed-Muller code 

RM(S 1 ,S 2 ,u,w 1 = 1,102 = si - u). (20) 
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If for an integer u' RM(S[, S 2 ,u' ,Wi,w' 2 ) is an optimal weighted Reed-Muller 



code of the same or smaller minimum distance as that of (20) then the latter 
code is of dimension at least that of the first one. 



Proof. If the latter code is of minimum distance close to that of (20) then 
it belongs to region II. As in the proof of the preceding proposition we have 
u' > us 2 j ' s' 2 . The improvement in dimension can be calculated to be at least 

/ ,wSl«2 S 2 U 
(S2-S 2 )(— — - -1) 

which takes on its minimal value s 2 /2 — 1 for u — s\ — 1. Combining this with 
the assumption s 2 > s 2 > 1 proves the proposition. □ 

Proposition 18. Consider positive integers si, s 2l s' 1; s' 2 with si-s 2 = s'iS 2) s[ > 
si > s 2 > s' 2 . Let u, s\ — 1 < u < (si — 1) + (s 2 — 1) be an integer and consider 
the optimal weighted Reed-Muller code 

RM{S 1 ,S 2 ,u,w 1 = l,w 2 = l). (21) 

If for an integer u' RM(S[, S' 2 ,u' ,w[ = l,w 2 ) is an optimal weighted Reed- 



Muller code of the same or smaller minimum distance as that of (21) then the 
latter code is of dimension at least that of the first one. 

Proof. The latter code either belongs to region II or III. For those in region 
III the result is pretty obvious so we consider only codes in region II. Let d be 



the minimum distance of the code in (21). We have v! > (sis 2 — d)/s 2 . The 



improvement in dimension can be calculated to be at least 

We may assume d > s' 2 as we are in region II and the result follows. □ 

The construction of weighted Reed-Muller codes is very concrete, but for 
completeness we should mention that it is not the most optimal. Consider 
instead the codes E(M,S) with S = Si x • • • x S m and 

M = {XI 1 ■ ■ ■ Xiz | ( Sl -»!)•■• (s m - i m ) > 6}. (22) 

Among the codes with designed distance 5 (Theorem |4| these are the codes of 
highest possible dimension. When S\ = ■ ■ ■ = S m = F q holds the construction 
simply is that of Massey-Costello-Justesen codes (see [^j and [TB]). 



4 Dual codes 

As is well-known, for the special case of S = F q x • • • x F q the duals of g-ary 
Reed-Muller codes, weighted Reed-Muller codes, and Massey-Costello-Justesen 
codes, respectively, are q-ary Reed-Muller codes, weighted Reed-Muller codes, 
and hyperbolic codes, respectively [28], [9] (for the definition of hyperbolic codes 
we refer to ( [24| below). More examples of codes E(M,S) where similar neat 
correspondences hold can be found in [3]. Turning to a general point ensemble 

5 = Si x • • • x S m , however, it does not in general hold that the dual of a 
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weighted Reed-Muller code is again a weighted Reed-Muller code. Nor does it 
hold in general that the dual of a Massey-Costello-Justesen code is a hyperbolic 
code. For a simple counter example which fits the description of a weighted 
Reed-Muller code as well as the description of a Massey-Costello-Justesen code 
consider the ordinary Reed-Solomon code over S = F* and recall that evs(l) is 
not a parity check for this particular code. 
Fortunately, for the class of codes 

E ± (M, S) ={ce Fg =lsl | c • ev s (M) = 0, for all M € M} 

we have a technique similar to that of Section [2] to estimate the minimum 
distance. This technique is known as the Feng- Rao bound ([7], [S]). We now 
recall this bound following the description of Shibuya and Sakaniwa in [27] , 
Consider the following definition of a linear code. 

Definition 19. Let B = {bx, • ■ • , b n } be a basis for and let G C B. We 

define C(B, G) = Span Fq {b \ b £ G}. The dual code is denoted C ± (B, G). 

The Feng-Rao bound calls for the following set of spaces. 
Definition 20. Let L_i = 0, L = {0} and L t = Span Fq {bi, . . . , &;} for I = 

1,.. .,71. 

We obviously have a chain of spaces {0} = L C L\ C • ■ ■ C L n _i C L n = 
F™. Hence, we can define a function as follows. 

Definition 21. Define p : F™ {0, 1, . . . , n} by p(v) — I if v e Lj\L/_i. 

Definition 22. Let I = {1,2,..., n}. An ordered pair € L 2 is said to be 
well-behaving if p(b u * b v ) < p{bi * bj) for all u and v with 1 < u < i,l < v < j 
and (u,v) ^ (i,j)- Here, * is the componentwise product. 

Definition 23. For I = 1, . . . , n define 

A*(0 = i^iihj) ^ I 2 I (hj) * s well-behaving and p(bi * bj) = 
The Feng- Rao bound now is ([571 Prop. 1]): 
Theorem 24. The minimum distance of C ± (B,G) is at least 

min{/2(0 | k i G}. 

Turning to the codes E(M,S) we enumerate the basis 

{evsiXl 1 ■■■X%) \ 0<ii< Sl ,...,0<i m <s m } 

according to a total degree lexicographic ordering on the monomials X] 1 ■ ■ ■ X^ . 
For bi = cvslXl 1 ■ ■ ■ X%p), < i\ < si,...,0 < i m < s m , we clearly have 
> • • • (*m+l)< Hence, by the Feng-Rao bound the minimum distance 

of E ± (M,S) is at least 

min{(ji + 1) • • • (i m + 1) | XI 1 ■ ■ ■ X 1 ™ £ M, < h < 8li . . . , < i m < s m }. (23) 
Consider the code 

(RM(<Si, . . . , S m , (s-l - 1) ■ ■ ■ (s m - I) - u - e,wx, . . . , Wm)) 1 - 
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where e is a very small positive number. The bound (23) tells us that the mini- 
mum distance is at least that of the weighted Reed-Muller code RM(5i, . . . , S m , 
u,w%, . . . ,w m ). Observe that the codes are of the same dimension. Similarly, 
the hyperbolic code E ± (M,S) 

M = {X^---X l ™\(i 1 + l)---(i rn + l)<5} (24) 

has designed minimum distance equal to S just as the Massey-Costello-Justesen 



code in (22). Again, the two codes are of the same dimension. 
The Feng-Rao bound comes with a decoding algorithm that corrects up to half 
the designed minimum distance [14]. This algorithm of course applies in 
particular to the above dual codes. 

The remaining part of the paper is concerned with decoding algorithms for the 
codes E(M,S) including the codes from Section [3] 



5 Subfield subcode decoding 



The Reed-Solomon code will be over the field F q ™ and is constructed by 

in the q m different elements 



As already noted by Kasami et al. in [17] , any ordinary q-ary Reed-Muller code 
(in the terminology of the present paper this means a g-ary Reed-Muller code 
from S = F q x ■ ■ • x F q ) can be seen as a subfield subcode of a Reed-Solomon 
code 

evaluating polynomials of degree at most uq m ~ 1 
of F q m. The above observation guarantees that codes E(M,S) in general can 
be seen as subcodes of subfield subcodes of certain Reed-Solomon codes over 
F q m , but it is not straightforward which elements of F q m to use. This problem, 
however, is easy to overcome if we use the approach by Santhi [25] . 
Let {&!,...,&„} be a basis for F 9 ™ as a vectorspace over F q . Following Santhi 

by 



we now define a map (p : ¥ q 



and note that <p(ai, . 
we have 



, a m ) = aibi H h a m b rn (25) 

= aib\ H \-a m b£. Writing X = cp(ai,. . . , a m ) 

bj] i 
b" 



b'i 





ai 




' x ' 








Xi 



















By [TH1 Cor. 2.38] the matrix on the left side is invertible and therefore there 
exist polynomials F\, . . . , F m g ¥ q m [T] such that ai = Fi(X). The polynomials 
Ft,..., F m do not depend on the point (ai, . . . , a m ) under consideration. From 
this observation one deduces that F(P) = F(F 1 ((p(P)), . . . , F m (ip(P))) for all 

p e f™. 

Theorem 25. Write S = {Pi,...,P n }. The code E(M,S) is a subcode of 
a subfield subcode of the Reed-Solomon code over F q m which is constructed by 
evaluating polynomials of degree at most 



t = max{deg M \ M e M} 



(26) 



in the elements <p(Pi 



,ip(P n ). Here ip is the function in (25). 
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Following the Pellikaan-Wu approach [3T] we can now decode E(M., S) by 
applying the Guruswami- Sudan list decoding algorithm to the corresponding 
Reed-Solomon code |12j and by performing a few additional steps. The com- 
plexity of the Guruswami-Sudan list decoding algorithm is in the literature often 
claimed to be C(n 3 ). For more precise statements of the decoding complexity 
which takes the multiplicity into account we refer to [3] . The Guruswami-Sudan 
algorithm corrects up to n(l — y/R) errors of the Reed-Solomon code. It is there- 
fore clear that the above approach can decode up to 



rn(l- / gm " 1 + 1 )l (27) 



errors of ^(M,^) (Here, t is as in (26)). This is indeed a fine result for many 



codes E(M,S). However, it is also clear that for other choices of _E(M,6>) (27) 



may be close to zero or even negative. For the particular case of an optimal 



weighted Reed-Muller code RM(Si, S2, u, wi — 1,102) (27) becomes 



r—^-V^ 1 - (28) 

(recall from Proposition [7] that w% < W2 always holds). If S1S2 is close to q 2 
and u is not too large then this bound guarantees list decoding. If s 2 is much 
smaller than si then the bound may not even guarantee that the algorithm can 
correct a single error. This is the reason why we in the present paper consider 
also a second decoding algorithm. Before getting to the second algorithm, we 
apply the first one to the Joyner code. 



5.1 The Joyner code 

Toric codes were introduced by Hansen in |13j and further generalized by Joyner 
in [15], by Ruano in (23J [24] and by Little et al. in [18]. Among the most famous 
toric codes is the [49, 11, 28] code over F 8 presented in [1^1 Ex. 3.9]. This code 
is known as the Joyner code. Attempts have been made to decode it, but with- 
out much luck so far. We now demonstrate how to decode it even beyond its 
minimum distance by applying the method of this section in combination with 
a small trick. 

The Joyner code originally was introduced in the language of polytopes. Alter- 
natively, one can define it [23] as a code E^M, S) where 

S = Fg x Fg = {(xi,yi), (2:49,2/49)}, 

M = {1} U {X l Y' J I 1 < i,j and i + j < 5}. 
Let the polynomial corresponding to a given code word c be 

F(X,Y) = F ,o+ F iJ xiYi - 

i,j > 1 
i + 3 < 5 

Let r — c + e be the received word. Assume for a moment that we know -Fb,o- 
We then subtract (-Fb.Oi ■ • • , ^0,0) from r to get a word r' that in the error free 
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Table 1: Error correction capability when the multiplicity used of the Reed- 
Solomon code decoder is m 



m 1 2 3 4 5 6 
capability 12 20 24 27 29 31 



positions corresponds to 

»,j>i 

i + i < 5 

For i = 1, ... ,49 we now divide the ith entry of r' with XiUi to produce a word 
7-". Observe, that this is doable because Xi,yi 7^ holds. The word r" in the 
error free positions corresponds to 

*,i > i 

i + j < 5 

We have r 7 ' = c" + e 7 ' where c" G E(M",S), M" = {X l Y^ \ i + j < 3} and e" 
is non-zero in exactly the same positions as e. The Reed-Muller code E(M",S) 
is a subfield subcode of a [49,25,25] Reed-Solomon code over F 64 . The exact 
form of the Reed-Solomon code is described by Theorem [25] Given a code word 
from the output of the Reed-Solomon list decoder we multiply for i = 1, . . . , 49 
the ith entry with x^yi and add -F^o- 

Of course we do not as assumed above know Fq^q in advance. Therefore we must 
try out all 8 possible values of this number. The error correction capability of 
the corresponding algorithm is described in Table [T] We see that we can correct 
up to 31 errors even though the minimum distance of the Joyner code is only 
28. 



6 An interpretation of the Guruswami- Sudan 
list decoding algorithm 

The second decoding algorithm of the present paper is a direct interpretation of 
the Guruswami-Sudan list decoding algorithm. We build on works by Pellikaan 
et al. [21], and Augot et al. [I], [2] who concentrated on g-ary Reed-Muller 
codes and Reed-Solomon product codes. We consider general code E(M 7 S) 
and improve on the above mentioned work by establishing new information on 
how many zeros of prescribed multiplicity a polynomial can have when given 
information about its leading monomial with respect to the lexicographic or- 
dering. In combination with a preparation step this will allow us to correct 
more errors. The idea of a preparation step comes from |10j . The improved 
information regarding the zeros is derived by strengthening results reported by 



Dvir et al. in [B]. This is done in Subsection 6.1 In Subsection 6.2 we present 



the algorithm and elaborate on its decoding radius. 
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6.1 Bounding the number of zeros of multiplicity r 

The definition of multiplicity that we will use relies on the Hasse derivative. 
Before recalling the definition of the Hasse derivative let us fix some notation. 
Assume we have a vector of variables X = (X 1} . . . ,X m ) and a vector k — 

(ki,...,k m )e Ng l then we will write X^ = X* 1 ■ ■ ■ Xf™ . In the following F is 
any field. 

Definition 26. Given F(X) £ F[X] and k £ N™ the k'th Hasse derivative of 
F, denoted by F^(X) is the coefficient of Z k in F(X + Z). In other words 

F(X + Z) = J2 F {i) {X)Z*. 

k 

The concept of multiplicity for univariate polynomials is generalized to mul- 
tivariate polynomials in the following way. 

Definition 27. For F(X) £ F[X}\{0} and a £ F m we define the multiplicity 
of F at a denoted by mult(F, a) as follows: Let M be an integer such that for 

every k = (fci, . . . , k m ) £ with k\ H h k m < M, F^(a) = holds, but 

for some k = (fci, . . . , k m ) £ with k\ H h k m = M, F&(a) ^ holds, 

then mult(F, a) — M . If F — then we define mult(F, a) = oo. 

The Schwartz-Zippel bound with multiplicity was reported already in pQ , [5] 
but was only recently proved, [S]. It goes as follows: 

Theorem 28. Let F(X) £ F[X] be a non-zero polynomial of total degree u. 
Then for any finite set SCF 

mult(F,a) < ulSl" 1 " 1 . 

We have the following useful corollary: 

Corollary 29. Let F(X) £ F[X] be a non-zero polynomial of total degree u 
and let S CF be finite. The number of zeros of F of multiplicity at least r from 
S m is at most 

-ISp- 1 . (29) 
r 

For the g-ary Reed-Muller codes 

RM q (u, m) = RM(F 9 , . . . , F q , u, 1, . . . , 1), 

Pcllikaan and Wu in [5T] presented two decoding algorithms, a subfield sub- 
code decoding algorithm and a direct interpretation of the Guruswami-Sudan 
algorithm. The analysis of the latter relies on [321 Lem. 2.4, Lem. 2.5] which 
combines to the following result: 

Proposition 30. Consider a polynomial F(X) £ F[X] of total degree u, u < rq 
and define w = [u/q\. The number of points in F™ where F has at least 
multiplicity r is at most equal to 

/m+r— 1\ rn i / \ /m+r— w — 2\ m—1 /m-\-r—w — l\ m 

{ m )q m + {u-qw)( m _j )q m - ( m )q m 



/m+r — 1\ 
V r-1 ) 
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Augot and Stepanov [T] gave an improved estimate on the decoding radius 
of the latter algorithm (the direct interpretation of the Gurswami-Sudan al- 



gorithm) by using instead Corollary 29 We here present a direct proof that 



indeed, Corollary [29] is stronger than Proposition [30| 

Proposition 31. For all u £ [0, rq — 1] it hold that \2!ty is smaller than or 
equal to (3C 



Proof. We consider the two expressions as functions in u on the interval [0,rg]. 



Our first observation is that ( 30 ) is a continuously piecewise linear function, 
each piece corresponding to a particular value of w. The corresponding r slopes 
constitute a decreasing sequence. Combining this observation with the fact that 



(291 is linear in u and with the fact that the two expressions are the same at 



each of the end points of the interval proves the result. □ 

As a preparation step to improve upon Theorem [28] and Corollary [29] we 
start by generalizing them. We will need a couple of results from [5] Sec. 2]. 
The first corresponds to Lem. 5]. 

Lemma 32. Consider F(X) £ F[X] and a £ F m . For any k = {k\, . . . , k m ) £ 



mult(F^\a) > mult(F, a) - (hi H h k m ). 

The next result that we recall corresponds to the last part of [6j Proposition 

6]. 

Proposition 33. Given F{X\, . . . , X m ) £ F[Xi, . . . , X m ] and 

Q(Y 1 ,...,Y l ) = (Q 1 {Y),...,Q m (Y))€F[Y 1 ,..., Y^ 1 

let F o Q be the polynomial F(Qi(Y), . . . , Q m (Y)). For any a £ F l we have 

mult(F o Q, a) > mult(F, Q(a)). 

We get the following Corollary, which is closely related to (S] Corollary 7]. 

Corollary 34. Let F(X 1 , . . . , X m ) £ F[X 1 ,...,X m ] and 6 X , . . . , 6 m _i, c € F m 
be given. Write F*(Ti, . . . ,T m _i) = F(T\bi + • • ■ + T! m _i& m _i + c). For any 
(fl, . . . , t m -i) £ F m_1 we have 

mult(F*{T u . . . ,T m _i), (ti, . . . ,t m -i)) 

> mult(F(Xi 7 . . .,X m ),txbi H htm-i&m-i + c). 

Let -< be the lexicographic ordering on the set of monomials in variables 
Xi, . . . , X m such that X m -< • • • -< X\ holds. We now write 

F(Xi,. . . , X m ) = X{ X ■ ■ ■ Xm-l Fh,-3m-i (Xm)- 

h, — ,jm-l 

Let XI 1 ■ ■ ■ XJ^ be the leading monomial of F with respect to -<. Then due to 
the definition of -<, Fi lt __i m l (X m ) is a (univariate) polynomial of degree i m . 
For a m £ F define 

r(a m ) = mvlt(F il! ... <im _ 1 (X m ),a m ). 
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Clearly, 

^ r(a m ) < i m . (31) 

a m es m 

We have 

jl, — ,3m-l 

and due to the definition of -< and to the definition of r{a m ) we have 

lm ^ (F (o,..,o,,(a m ) )(Xij _ > x m _ x , a m )) = X\- ■ ■ ■ X 1 ™Z{ . (32) 



Applying first Lemma 32 with k = (0, . . . , 0, r(a m )) and afterwards Corollary 34 
with bi = (1, 0, . . . , 0), . . . , 6 m -i = (0, . . . , 0, 1, 0), c = (0, . . . , 0, o m ) and t x = 
a%, . . . ,t m — l — a m-i we get the following result which is closely related to a 
result in [6j Proof of Lemma 8] : 

mult(F(Xi, . . .,X m ), (ai, . . .,a m )) 

< (0 + • • • + + r(a rn )) + mult(F(°-- °' r ^(X lt X m ), (a u . . . , a m )) 

< r{a m ) + mvlt(F^-'°' r ^\X u . . .,X m _ x ,a m ), (a u . . . , a m _ x )). (33) 

We are now ready to generalize Theorem |28| Let in the remaining part of this 
subsection Si, ... , S m be finite subsets of arbitrary field F. Also we will relax 
from the assumption that si > ■ ■ • > s rn . 

Theorem 35. Let F(X) E F[X] be a non-zero polynomial and let lm(F) — 
XY ■ ■ ■ X % ™ be its leading monomial with respect to a lexicographic ordering. 
Then for any finite sets S\, ... , S m C F 

^ mult(F, a) < i\s 2 ■ ■ ■ s m + sii 2 S3 • ■ ■ s m H h s% ■ ■ ■ s m ^ x i m . 

a£Si X ■■■XS m 

Proof. We prove the theorem for the monomial ordering -<. Dealing with gen- 
eral lexicographic orderings is simply a question of relabeling the variables. 



Clearly the theorem holds for m = 1. For m > 1 we consider (33). Assuming 
the theorem holds when the number of variables is smaller than m we get by 
applying |3T|) and (p2l the following estimate 



mult(i r , a) 

< i m si ■ ■ ■ s m -i + s m (iis 2 ■ ■ ■ s m -i H h i m _isi • • • s m _ 2 ) 

= hs2 ■ ■ ■ s m + i2Sis 3 ■ ■ ■ s m H i m si ■ ■ ■ s m ~i 

as required. □ 

We have the following immediate generalization of Corollary [29} 

Corollary 36. Let F(X) 6 F[X] be a non-zero polynomial and let lm(F) — 
X^ 1 ■ ■ ■ X 1 ™ be its leading monomial with respect to a lexicographic ordering. 
Assume S%, ... , S m C F are finite sets. Then over Si x • • • x S m the number of 
zeros of multiplicity at least r is less than or equal to 

(hs 2 ■■■s m + sii 2 s 3 ■ • • s m H h sx ■ ■ ■ s TO _i« ro )/r. (34) 
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The analysis leading to Theorem [28] suggests the following function to more 
accurately estimate the number of zeros of multiplicity at most r of a polynomial 
with leading monomial X\ x ■ ■ ■ X^ 1 : 

Definition 37. Let r € N,«i, . . . , i m G N . Define 

D(ii,r,si) = min{[^-J,si} 

and for m > 2 

D{i\, . . . , i m , r, Si, . . . , s m ) 



max < (s m - m u r )D(ix, . . ,,i m -i,r,si, . . . ,s m -i) 

(ui ; ...,u r )£A(i m ,r,s m ) ^ 

+ uiD(ii, . . . ,i m -x,r - 1, si, . . . , s m _i) H 

+ U r -1-D(il, . . . , i m _i, 1, Si, ... , s m _i) + u r s x • • • S m _! 

{(wi, . . . , u r ) € Nq I u\ H — • + u r < s m a77,rf Mi + 2it2 + • • • + ru r < z m }. 

Theorem 38. For a polynomial F(X) £ F[X] let Xf 1 •• • X^ be its leading 
monomial with respect to -< ( this is the lexicographic ordering with X m -<•••-< 
Xij. TTien F has at most D(ii,...,i m ,r,si,...,s m ) zeros of multiplicity at 
least r in Si x ■ ■ ■ x S m . The corresponding recursive algorithm produces a 
number that is at most equal to the number found in Corollary \3b] and is at 
most equal to s\ ■ ■ ■ s m . 

Proof. The proof of the first part of the proposition is an induction proof. The 
result clearly holds for m = 1. Given m > 1 assume it holds for m — 1. For 
d = 1, ... ,r — 1 let Ud be the number of a m 's with r(a m ) = d and let u r be 
the number of a m 's with r(a m ) > r. The number of a m 's with r(a m ) = 
is s m — u\ — ■ ■ ■ — u r . The boundary conditions that u\ + • • • + u r < s m 
and u\ + 2i*2 + • • • + ru r < i rn are obvious. For every a m with r(a m ) = d, 
d = 0, . . . , r — 1 for (ai, . . . , a m ) to be a zero of multiplicity at least r the 



last expression in (33) must be at least r — d. For a m with r(a m ) > r all 
choices of oi, . . . , a m -i are legal. This proves the first part of the proposition. 
As both Corollary 36 and the above proof rely on (33), Theorem 38 cannot 
produce a number greater than what is found in Corollary |36| The condition 
Ui + ■ ■ ■ + u r < s m and the definition of D(i\, r, s\) imply the last result. □ 

It only makes sense to apply the function D(ii, . . . , i m , r, s±, . . . , s m ) to 
monomials X] 1 ■ ■ ■ X^ in 

A(r, si, . . . , a m ) = {XI 1 ■ ■ ■ X^ | [ii/sij + • • • + bm/smj < r}. 

Proposition 39. Assume X'^ 1 ■ ■ ■ X]™ $ A(r, s l5 . . . , s m ). Then there exists a 
polynomial with leading monomial X 1 ^ ■ ■ ■ X*™ such that all elements of Si x 
■ ■ ■ x S m are zeros of multiplicity at least r. 
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Example 40. In a number of experiments listed in \ 11]/ we calculated the value 
D(ii, . . . , i m , r, q, . . . , q) for various choices of m, q and r and for all values 
of . . . , i m ) such that X\ m ■ ■ ■ X^ 1 € A(r, q, . . . , q). Here we list the mean 
improvement in comparison with the situation where Corollary \2S\ is applied. 
More formally, we list in Table^for various fixed q,r,m the mean value of 

min{(ii j h i m )g m ~ 1 /?-, q m } - D(h, . . . ,i m ,r,q, . . . ,q) 

mm{(i 1 + ---+i m )q m - 1 /r,q m } 

Despite the significant mean improvement, according to our experiments in 111}/ 



Table 2: The mean value of (351; truncated 



0.363 
0.217 
0.191 
0.155 
0.128 
0.126 



0.273 
0.286 
0.197 
0.167 
0.137 
0.127 



0.337 
0.228 
0.232 
0.174 
0.138 
0.134 



0.291 
0.236 
0.195 
0.197 
0.138 
0.126 



0.301 
0.194 
0.158 
0.139 
0.119 
0.114 



0.300 
0.224 
0.169 
0.145 
0.122 
0.115 



0.342 
0.213 
0.180 
0.148 
0.121 
0.113 



0.307 
0.214 
0.172 
0.153 
0.119 
0.111 



0.248 
0.158 
0.125 
0.110 
0.093 
0.089 



0.260 
0.177 
0.135 
0.116 
0.098 
0.093 



for most fixed degrees u there are examples of exponents . . . , i m ), i\ + • • ■ + 
im = u such that D(ii, . . . ,i m ,r, si, . . . , s m ) = H h i m )q m ~ 1 /r\ . 

Sometimes the values -D(«i, . . . , i m , f, Si, . . . , s m ) may be time consuming to 
calculate. Therefore it is relevant to have some closed formula estimates of these 
numbers. We next present such estimates for the case of two variables. Note, 
that the following proposition covers all monomials in A(r, s\, S2). 

Proposition 41. For k = 1, . . . ,r — X, D(i±, 12, r, sx, S2) is upper bounded by 
(C.1J s 2 h 



r r r — k 

if (r - fc)^jSi < i\ < (r — k)si and < i 2 < ks 2 , 

(C.2) s 2 ^ + ((k + 1) 82 - i 2 ){^j --) + (i a - ks 2 )( Sl - ?i) 
r r — k r r 

if (r — k)^rp[Si < %\ < (r — k)si and ks 2 < i 2 < (k + l)s 2 , 

(C3) ^+ 

r k + 1 r 
if (r — k — l)s% < i\ < (r — fc)^rjSi and < i 2 < (k + l)s 2 - 

Finally, 

(C.4) D(i 1 ,i 2 ,r,s 1 ,s 2 ) = s 2 [—\ +i 2 (si - [— J) 

r r 
if Si(r — 1) < i\ < sir and < i 2 < s 2 . 

The above numbers are at most equal to min{(iiS2 + s\i 2 )/r, sis 2 }. 

Proof. First we consider the values of i\,i 2 , r, si, s 2 corresponding to one of the 
cases (C.l), (C.2), (C.3). Let k be the largest number (as in Proposition 41) 
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such that i\ < (r — k)s\. Indeed k € {1, . . . , r — 1}. We have 



D(i 1 ,i 2 ,r, s 1 ,s 2 ) < 

( k . , k k, . . , k 

max <s 2 \-Ui{ )H \- u k ( ) 

(u 1 ,...,u r )£B{i 2 ,r,s 2 ) [ r r — 1 r r—kr 

+ u fe+ i(si ) H htt r (si ) \ (36) 

r r J 

where 

B(i 2 ,r,s 2 ) = {(ux,...,u r ) € Q r | < u lt . . . ,u r , 

itl H \-u r < s 2 ,u 1 + 2u 2 + • ■ ■ + ru r < i 2 }. 

We observe, that 

*(-V ^) < ^) 

r — t r r — k r 

holds for I < fc. Furthermore, we have the biimplication 

(r - fc)-^! < k (fc + - -) > - -)• 

r + 1 r — kr r 



Therefore, if the conditions in (C.l) are satisfied then (36) takes on its maximum 



when Uk = ¥■ and the remaining Uj's equal 0. If the conditions in (C.2) are 



satisfied then (36) takes on its maximum at u k — {k + l)s 2 — i 2l iik+i = (i 2 — ks 2 ) 



and the remaining itj's equal 0. If the conditions in (C.3) are satisfied then (36) 
takes on its maximal value at Uk+i = -j^j and the remaining itj's equal 0. 
Finally, if si(r — 1) < k < s i r an d < 12 < s 2 then D(k, ki r i s i> s 2) is the 
maximal value of 

s 2 L-J+«i(si-L-J) + --- + ^(.<q-L-J) 

r r r 

over -B(i2, r, s 2 ). The maximum is attained for u\ = i 2 and all other itj's equal 0. 
The proof of the last result follows the proof of the last part of Theorem [38] □ 

Remark 42. Experiments show (see 111]) that the numbers produced by Propo- 



sition 41 are often much smaller than rrrm{(iiS2+si^2)A*) sis 2 }. However, there 
are cases where they are identical. This happens for example when i\ = Si(r— 1) 
and r divides s 1 ands 2 . In the proof of (C.l), (C.2), (C.3) we allowed u±, ... ,u r 
to be rational numbers rather than integers. Therefore we cannot expect the 



upper bounds in Proposition 41 to equal the true value of D{ii,i 2 ,r, s\, s 2 ) in 
general. Our experiments show that the bounds in (C.l), (C.2), (C.3) are some- 
times close to D(ii,i 2 ,r, s\, s 2 ) but not always. Hence the best information is 
found by actually applying the function D{i\,i 2 ,r, S\, s 2 ) directly. 

6.2 The decoding algorithm 

The main ingredient of the decoding algorithm is to find an interpolation poly- 
nomial 

Q(Xi, . . . , X m , Z) = Qq{X\, . . . , X m )+Qi(Xi, . . . , X m )Z+- ■ -+Q t (Xi, . . . , X^Z 1 
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such that Q{X\, . . . , X m , F{X\, . . . , X m )) cannot have more than n—E different 
zeros of multiplicity at least r whenever Supp(F) C M. The integer E above 
is the number of errors to be corrected by our list decoding algorithm. In [21] , 
PQ , [2] this requirement is described in terms of bounds on the total degree of 
the polynomials Qj. As we will use improved information that depends not 
on total degree but on the leading monomial with respect to a lexicographic 
ordering the situation becomes more complicated. To fulfill the above require- 
ment we will define appropriate sets of monomials B(i, E,r), i — l,...,t and 
then require Qi(Xi, , . . ,X m ) to be chosen such that Supp(Qi) C B(i,E,r). 
Rather than using the results from the previous section on all possible choices 
of F{X\, . . . ,X m ) with Supp(i^) C M we need only consider the worst cases 
where the leading monomial of F is contained in the following set: 

Definition 43. 

M = {M e M | if N e M and M\N then M = N}. 

Hence, M is so to speak the border of M. 

Definition 44. Given positive integers i,E,r with E < n let 

B(i, E, r)={K G A(r, s u . . . , s m ) \ D r {KM l ) <n- E for all M G M}. 

Here D r (X' L 1 , . . . , X%?) can either be D(i\, . . . , i m , r, s\, . . . , s m ) or in the case 
of two variables it can be the numbers from Propo siti on \41\ Another option 



would be to let D r {X^ , . . . ,X^) be the number in {34) 



The decoding algorithm calls for positive integers t, E, r such that 
t 

J2\B(i,E,r)\>nN(m,r), (37) 

where N(m, r) — (™^) is the number of linear equations to be satisfied for a 
point in F™ +1 to be a zero of Q(Xi, . . . , X m , Z) of multiplicity at least r. As 



we will see condition (37 1 ensures that we can correct E errors. We say that 



(t,E,r) satisfies the initial condition if given the pair (E,r), t is the smallest 



integer such that (37) is satisfied. Whenever this is the case we define B'{t, E, r) 



to be any subset of B(t, E, r) such that 



t-i 

\B(i,E,r)\ + \B'(t,E,r)\ = nN(m,r) + 1. 

i=i 

Replacing B(t, E, r) with B'(t, E, r) will lower the run-time of the algorithm. 

Algorithm 1. Input: 

Received word r = (r\, . . . , r n ) £ F™. 

Set of integers (t, E, r) that satisfies the initial condition. 

Corresponding sets B(l, E,r) . . . , B(t — 1, E, r), B'(t, E, r). 

Step 1 

Find non-zero polynomial 

Q(X\, . . . , X m Z) — Qo(Ai, . . . , X m )+Qi(Xi, . . . , X m )Z+- ■ -+Q t (Xi, . . . , X m )Z l 
such that 
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1. Supp{Qi) C B(i, E,r)fori = l,...,t-l and Supp(Q t ) C £, r), 

2. (Pj,rj) is a zero of Q(X\, . . . , X m , Z) of multiplicity at least r for i = 
l,...,n. 

-find a/? F(Xi, . . . , X m ) £ ¥ q [X\, . . . , X m ] such that 

(Z - F(X U . . . ,X m ))\Q{X u . . .,X m , Z). (38) 

Output: 

A list containing (F(P±), . . . ,F(P n )) for all F satisfying (38). 



Theorem 45. The output of Algorithm^ contains all words in E(M,S) within 
distance E from the received word r. Once the preparation step has been per- 
formed the algorithm runs in time 0(n 3 ) where n = n(™^). ^or given mul- 
tiplicity r the maximal number of correctable errors E and the corresponding sets 
5(1, E,r),..., B(t-1, E, r), B'(t, E, r) can be found in time 0(n log(n)r"V|M|/cr) 
assuming that the values of the function D r are known. Here a = maxjdeg M \ 
M G M} and s' = max{si, . . . , s m }. 

Proof. The interpolation problem corresponds to n homogeneous linear equa- 
tions in n + 1 unknowns. Hence, indeed a suitable Q can be found in time 
C(n 3 ). Now assume Supp(P) C M and that dist H (ev s (F),r) < E. Then 
Pj is a zero of Q(X\, . . . , X m , F(Xi, . . . , X m )) of multiplicity at least r for 
at least n — E choices of j. By the definition of B(i,E,r) this can, how- 
ever, only be the case if Q{X\, . . . , X m ,F(Xi, . . . , X m )) = 0. Therefore, Z — 
F(Xi, . . . , X m ) is a factor in Q(Xi, . . . , X m , Z). Finding linear factors of poly- 
nomials in (¥ q [X\, . . . ,Xm\)[Z] can be done in time C(n 3 ) by applying Wu's 
algorithm in [2S] (see [22} p. 20]), □ 

Algorithm[l]works for general codes E(M, S) and for any of the three possible 
choices of D r (Xl 1 ■ ■ ■ X]^ ) as described prior to the algorithm. In such a general 
setting it is impossible to say anything reasonable regarding the decoding radius. 
The algorithm apparently works best for not too large code dimensions. With 
this in mind we restrict the analysis to optimal weighted Reed-Muller codes 
RM(jSi, 52, u, w\, W2) in region I. That is, we assume w± = 1, W2 = S1/S2 
and u < Si — Si/s2- As the function D(ii, 12, r, Si, S2) is highly irregular and 
Proposition [41] contains four quite different cases it seems impossible to perform 
the analysis for other choices than D{i\, 12) = {i\S2 + «2Si)/r which corresponds 
to the weakest version of the decoding algorithm. 

Proposition 46. Consider an optimal weighted Reed-Muller code RM{S\ 1 S2, u, ui\ 
1, W2 = S1/S2) with S2\si and u < S\ ~ S1/S2 a positive integer. When equipped 
with D(ii,i2) — {i\S2 + i2S\)/r the decoding radius of Algorithm^ is at least 

sis 2 (l- (39) 
Proof. Let v be divisible by u. The number of variables in the interpolation 
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polynomial when t = deg z Q is chosen to be v/u is lower bounded by 



v/u—l 

£ 

3=0 



/u-1 

£ [( V + 1 ~ 3 U ) + o( U + 1 ~ ^ J«) s 2/Sl - 1)] 



1 ' u / u 3 

> ~—y{uif>—~. 

2 si si 6u 

2 — 1 

The number of equations is Sis 2 r(r + l)(r + 2)/6 and therefore 



u > \Jur{r + l)(r + 2)s\ 
is a sufficient condition for the existence of an interpolation polynomial. Assume 
E < s lS2 (l - ^/u(l + l/r)(l + 2/r)/si) 

E < sis 2 - ^/w(r + l)(r + 2)sf. 

Substituting u = \J ur(r + f )(r + 2)s^ we get r(siS2 — E) > vs 2 which ensures 
that Q(Xi, X2(F(Xi, X2)) = for any codeword c — cvs(F) within distance E 
from r. Letting r go to infinity finishes the proof. □ 



Comparing the decoding radii ( 28 1 and ( 39 ) we conclude that when s 2 is 



close to q then the subfield subcode decoder is superior. On the other hand 
when S2 is much smaller than q then the decoding algorithm of the present 
section performs best. 

Example 47. In this example we investigate the performance of Algorithm [7] 
when applied to optimal weighted Reed-Muller codes and Massey-Costello-Justesen 
codes coming from the point ensembles S = S\ X S2 with s± = 64, s 2 = 
8, and s± — 256, s 2 = f 6, respectively. Our findings are presented in Ta- 
ble [3| and Table [^J respectively. The decoding capability is calculated for dif- 
ferent choices of D^X^ 1 X l 2 2 ) and different multiplicities r. The symbol S, C, 
and D, respectively, corresponds to D r (X' L 1 X^ 2 ) being chosen as the Schwartz- 



Zippel bound \3Jj\ , the closed formulas of Proposition 4F and the function 
D(i2,ii,r,S2,si), respectively. The letter W stands for optimal weighted Reed- 
Muller code and I means the Massey-Costello-Justesen code of the same mini- 
mum distance. Further u is the third argument in the notion RM(S\, S2, u, W\ = 
1,102 = si/s 2 ) and d is the minimum distance. Sub stands for the estimated 



decoding radius (21) of the algorithm in Section [5| and Dim is the dimension 



of the code. For large values of r the calculations regarding D(ix, 12, r, s%, S2) 
become quite heavy and have therefore not been made. We can see from the 
tables that for the considered codes Algorithm^ outperforms the subfield subcode 
approach from Section [3[ In some cases it decodes much more than half the 
minimum distance. It is apparent that the function D{ii, i 2 , r, s±, s 2 ) as well 
as the closed formula expressions of Proposition help bringing up the error 
correction capability in comparison with the situation where the Schwartz-Zippel 



bound {34) is used. It is clear that the small gain in dimension by considering 



Massey- Costello-Justesen codes rather than optimal weighted Reed-Muller codes 
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Table 3: Table of error correction capability for optimal weighted Reed-Muller 
codes and Massey-Costello-Justesen codes when s\ = 64 and = 8. 





ii Id 

LL j LL 






7 4^fi 


1 u 


309 


1U 




90 


"^9 


r 


Uound 


W I 


W I 


W I 


W 


I 


W 


I 


w 


I 




S 


267 


243 


191 


103 


95 


95 


87 


67 


59 


i 
A 


? < 
L 


2»6 


26b 


1 1 n 
219 


131 


1 no 

12o 


122 


iiy 




94 




D 


298 


277 


228 


135 


131 


121 


119 


99 


95 





S 


287 


263 


213 


130 


122 


122 


117 


95 


90 


3 


c 


301 


279 


234 


149 


145 


138 


135 


113 


109 




D 


319 


298 


255 


177 


175 


161 


160 


139 


135 




S 


295 


273 


225 


145 


139 


139 


131 


111 


105 


4 


c 


307 


286 


242 


159 


155 


147 


145 


123 


118 




D 


328 


311 


269 


196 


195 


181 


181 


160 


159 




s 


312 


292 


247 


173 


166 


166 


159 


140 


134 


9 


c 


318 


299 


255 


178 


173 


169 


166 


144 


139 




s 


320 


301 


258 


185 


178 


178 


171 


153 


147 


20 


c 


323 


304 


262 


188 


182 


180 


175 


155 


149 




Sub 


198 


149 


33 





















243 


239 


227 


195 


191 


175 




Dim 


4 


5 


8 


24 


25 


27 


28 


39 


41 



comes with a heavy price as Algorithm [7] corrects much fewer errors. By in- 
spection the estimation of decoding radius from Proposition \46\ seems to be quite 
close to what is found by our computer experiments. 

7 Conclusion remarks 

In this paper we have shown that weighted Reed-Muller codes are much better 
than their reputation when defined over general point ensembles S = S± x • • • x 
S m . We treated in detail the case m — 2 and gave some results for m > 2. It 
is a subject of future studies to also establish detailed information for the case 
m > 2. We derived two decoding algorithms that work well for different classes 
of weighted Reed-Muller codes and affine variety codes E^M, S) in general. For 
not too high dimensions these algorithms perform list decoding. For higher 
dimensions it is a subject of future research to design list decoding algorithms. 
Using the first algorithm in combination with some extra operations we decoded 
the [49, 11,28] Joyner code beyond its minimum distance. It is apparent that 
such an approach would work for other toric codes coming from polytopes of 
the same shape. 

This work was supported in part by Danish Natural Research Council grant 272- 
07-0266. The authors gratefully acknowledge support from the Danish National 
Research Foundation and the National Natural Science Foundation of China 
(Grant No. xxx) for the Danish-Chinese Center for Applications of Algebraic 
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Table 4: Table of error correction capability for optimal weighted Reed-Muller 
codes and Massey-Costello-Justesen codes when si — 256 and s 2 = 16. 





u/d 


5 4016 


8 3968 


15 3856 


31 


3600 


36 


3620 


55 


3216 


r 


Bound 


W I 


W I 


W I 


W 


I 


W 


I 


W 


I 


— 


S 


2591 


2335 


1927 


1359 


1335 


1231 


1207 


839 


791 


2 


c 


2680 


2456 


2112 


1565 


1557 


1392 


1391 


1022 


1003 




D 


2729 


2504 


2153 


1589 


1583 


1411 


1408 


1035 


1015 




S 


2714 


2479 


2106 


1578 


1551 


1455 


1434 


1082 


1034 


Q 
O 


c 






9940 


1695 


1684 


1552 


1547 


1190 


1167 




D 


2861 


2651 


2326 


1859 


1855 


1707 


1706 


1359 


1351 




S 


2779 


2555 


2195 


1691 


1667 


1575 


1551 


1211 


1163 


4 


c 


2843 


2635 


2305 


1782 


1767 


1638 


1632 


1284 


1260 




s 


2894 


2689 


2362 


1895 


1871 


1784 


1763 


1443 


1367 


9 


c 


2928 


2730 


2415 


1935 


1919 


1811 


1804 


1469 


1442 




s 


2947 


2751 


2439 


1988 


1966 


1882 


1862 


1551 


1506 


20 


c 


2964 


2772 


2464 


2007 


1989 


1894 


1884 


1562 


1529 




Sub 


1806 


1199 


130 





















2007 


1983 


1927 


1799 


1759 


1607 




Dim 


4 


5 


8 


24 


25 


27 


28 


39 


41 
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